A Study of D° -> K° S K° S X Decay Channels 



The FOCUS Collaboration* 

J. M. Link a , P. M. Yager a , J. C. Anjos b , I. Bediaga b , 
C. G6bel b , A. A. Machado b , J. Magnin b , A. Massafferri b , 
J. M. de Miranda 13 , I. M. Pepe b , E. Polycarpo b , 
CN ! A. C. dos Reis b , S. Carrillo c , E. Casimiro c , E. Cuautle c , 

^! A. Sanchez-Hernandez c , C. Uribe c , F. Vazquez , L. Agostino d , 

L. Cinquini d , J. P. Cumalat d , B. O'Reilly d , I. Segoni d , 



O 

(N 



K. Stenson d , J. N. Butler 6 , H. W. K. Cheung 6 , G. Chiodini 6 , 
I. Gaines 6 , P. H. Garbincius 6 , L. A. Garren 6 , E. Gottschalk 6 , 

P. H. Kasper 6 , A. E. Kreymer 6 , R. Kutschke 6 , M. Wang 6 , 
L. Benussi f , M. Bertani f , S. Bianco f , F. L. Fabbri f , A. Zallo f , 
M. Reyes g , C. Cawlfield h , D. Y. Kim h , A. Rahimi h , J. Wiss h , 
2 ! R- Gardner 1 , A. Kryemadhi 1 , Y. S. Chung j, J. S. Rang-", 

B. R. Kc-j, J. W. KwakJ, K. B. Lee j , K. Cho k , H. Park k , 
G. Alimonti^, S. Barberis^, M. Boschini^, A. Cerutti^, 
P. D'Angek/, M. DiCorato^, P. Dim £ , L. Edera^, S. Erba £ , 
P. Inzani^, F. Leveraro^, S. Malvezzi^, D. Menasce^, 
M. Mezzadri^, L. Moroni^, D. Pedrini^, C. Pontoglio^, 
F. Prelz £ , M. Rovere^, S. Sala £ , T. F. Davenport IIP 1 , 



X 



V. Arena 11 , G. Boca 11 , G. Bonomi n , G. Gianini 11 , G. Liguori 11 , 
D. Lopes Pegna 11 , M. M. Merlo 11 , D. Pantea 11 , S. P. Ratti 11 , 
C. Riccardi 11 , P. Vitulo n , H. Hernandez , A. M. Lopez , 
H. Mendez , A. Paris , J. Quinones , J. E. Ramirez , 
Y. Zhang , J. R. Wilson p , T. Handler q , R. Mitchell q , 
D. Engh r , M. Hosack 1 ', W. E. Johns r , E. Luiggi r , J. E. Moore r , 
M. Nehring r , P. D. Sheldon 1 ", E. W. Vaandering r , M. Webster 1 ", 

M. SheafF 

a University of California, Davis, CA 95616 
h Centro Brasileiro de Pesquisas Fisicas, Rio de Janeiro, RJ, Brasil 
C CINVESTAV, 07000 Mexico City, DF, Mexico 



See http://www-focus.fnal.gov/authors.html for additional author information. 



Preprint submitted to Elsevier Science 



7 February 2008 



d University of Colorado, Boulder, CO 80309 
e Fermi National Accelerator Laboratory, Batavia, IL 60510 
f Laboratori Nazionali di Frascati dell'INFN, Frascati, Italy 1-00044 
g University of Guanajuato, 37150 Leon, Guanajuato, Mexico 
^University of Illinois, Urbana- Champaign, IL 61801 
^Indiana University, Bloomington, IN 47405 
J Korea University, Seoul, Korea 136-701 
^Kyungpook National University, Taegu, Korea 702-701 
^INFN and University of Milano, Milano, Italy 
m University of North Carolina, Asheville, NC 28804 
n Dipartimento di Fisica Nucleare e Teorica and INFN, Pavia, Italy 
° University of Puerto Rico, Mayaguez, PR 00681 
p University of South Carolina, Columbia, SC 29208 
q University of Tennessee, Knoxville, TN 37996 
r Vanderbilt University, Nashville, TN 37235 
s University of Wisconsin, Madison, WI 53706 



Abstract 



Using data from the FOCUS experiment (FNAL-E831), we report on the decay 
of D° mesons into final states containing more than one Kg. We present evidence 
for two Cabibbo favored decay modes, D° -» K%K%K~ir + and D° -» K%K%K + ir~, 
and measure their combined branching fraction relative to D° — > K®it + it~ to be 

0.0106 ± 0.0019 ± 0.0010. Further, we report new measure- 



r(D°-^/t°/t° ^ ± 7rT) 



r(D°^^°7r+7r-) 

ments of ^Z^fj^ = 0-0179 ± 0.0027 ± 0.0026, JlgggL = .0144 ± 

0.0032 ± 0.0016, and ^^illl-^ = °- 0208 ± - 0035 ± °- 0021 where the first 
error is statistical and the second is systematic. 



1 Introduction 



Detailed measurements of rare exclusive decay modes of charm mesons provide 
a powerful way to probe the details of charm decay such as the contributions 
of W exchange diagrams and final state interactions. These measurements aid 
our understanding of the interplay between the weak and strong interactions, 
mainly for multibody decays where theoretical predictions are poorer than for 
two body decays. We report the first observation of D° — > KgKgir ± K T , and 
new measurements of D° -> K° S K° S K° S , D° -> K° S K° S , and D° -> K° s K° s ir + 7i- . 
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The fixed-target charm photoproduction experiment FOCUS, an upgraded 
version of E687 [1], collected data during the 1996-1997 fixed-target run at 
Fermilab. A photon beam is derived from the bremsstrahlung of secondary 
electrons and positrons produced from the 800 GeV/c Tevatron proton beam. 
The photon beam interacts with a segmented beryllium-oxide target. The av- 
erage photon energy for reconstructed charm events is 180 GeV. 
Two silicon microvertex systems provide excellent separation between the pro- 
duction and charm decay vertices. One silicon strip detector called target 
silicon (TS) is embedded in the BeO target segments [2]. The other silicon 
strip detector (SSD) is located downstream of the target region. Charged par- 
ticles are tracked and momentum analyzed with five stations of multiwire 
proportional chambers in a two magnet forward spectrometer. Three multi- 
cell threshold Cerenkov detectors are implemented to identify electrons, pions, 
kaons, and protons [3]. 



2 Reconstruction of individual particles: K s , K, and 7r 



A detailed description of individual K s reconstruction in the FOCUS spec- 
trometer can be found elsewhere [4]. Briefly, Kg's are identified in several 
different regions depending on the K s decay length; upstream of the magnet 
with silicon strip information, upstream of the magnet without silicon strip 
information, and inside of the magnet. Depending on whether the daughter 
pions from the K s pass through the second magnet, and hence have a well de- 
fined momentum, a series of different techniques are employed. Roughly 15% 
of the K s are found upstream of the silicon system and can be used in locat- 
ing the D° decay vertex. When there are two or more K s in a decay channel, 
the daughter tracks must all be distinct and can not be shared between K s 
candidates. The reconstructed mass of the K s must be within three standard 
deviations of the nominal K s mass. 

All charged tracks from the charm decay must be singly linked to the sili- 
con microstrip, be of good quality, and inconsistent with zero degree photon 
conversion. Cerenkov particle identification (PID) for charged particles is per- 
formed by constructing a log likelihood value Wj for the particle hypotheses 
(i = e, 7r, K, p). The ir consistency of a track is defined by AW W = W mi „ — 
where W m in is the minimum W value of the other three hypotheses. Similarly, 
we define A Wk,k = — W/r and AWk,p = VV P — Wk for kaon identification 
and we require AWk,tt > 1 and AWk, p > —2 for kaons and AW n > —5 for 
pions. 
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3 Reconstruction of D° candidates 



The D° candidates are formed by making a vertex hypothesis for the daugh- 
ter particles whenever possible. We use the K s candidates and two charged 
tracks of the correct charge combination to form a D° candidate. The confi- 
dence level of the decay vertex of the D° candidate is required to be greater 
than 1%. The combined momentum vector located at the decay vertex forms 
the D° track. Using the D° track as a seed track for the candidate driven ver- 
texing algorithm [1], we search for a production vertex in which the confidence 
level is greater than 1% and the primary multiplicity including the seed track 
must have at least 3 tracks. The production vertex must be inside the target. 
The significance of separation between the production and the decay vertices 
(L/a L ) must be greater than 4 for D° — > K s K s tc t K ± and greater than 6 for 
D° — > KgKgir + iT~ . Different values are chosen for the L/a L cut due to the 
level of background for the two channels, but scans are performed in L/ol as 
part of the systematic uncertainty studies. 

A stand-alone vertex algorithm is used to reconstruct the primary vertex of the 
event for decay modes without charged tracks [1]. All the silicon microvertex 
tracks in the event, excluding those already assigned to pions used to recon- 
struct Kg, are used to form all possible vertices of the event with a confidence 
level greater than 1%: we then choose the primary vertex to be the highest 
multiplicity vertex. Ties are resolved choosing the most upstream vertex. 

The signal channels are normalized to D° — > K 7r + 7r~ which is the most 
abundant D° decay mode containing a K s . We use only the most basic cuts 
in this analysis. The invariant mass distribution for D° — > Alj7r + 7r~ at an 
L/aL > 4 is presented in Fig. la. To determine the reconstruction efficiency of 
this channel we use a Monte Carlo simulation using our best determined de- 
cay distribution to generate D° K°n + n~ . The reconstructed Monte Carlo 
sample was six times larger than the data sample. We fit our signals from both 
data and Monte Carlo with a Gaussian for the signal and a 2 nd degree polyno- 
mial for the background using a maximum likelihood fit. As our most copious 
K s decay channel, we use the D° — > K s ir + ir~ mode to search for variations 
versus reconstructed K s categories. We identify a systematic uncertainty on 
K s reconstruction of 7.1%. 

The systematic uncertainty for each branching fraction is independently de- 
termined. Several tests are performed. To check the stability of the branching 
fractions we fix the cuts to their standard values and vary one cut at a time. 
Cuts which are scanned include particle identification, significance of separa- 
tion, isolation of the secondary vertex from any other track, and confidence 
level of the secondary vertex. We also perform split sample systematics by 
dividing our samples into two sub-samples based on momentum and data- 
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taking period. We study fit variant systematics by calculating the branching 
fraction when the signals are fit with 1 st , 2 nd , or 3 rd degree polynomials. 

We look for Monte Carlo resonance systematics by searching for a difference 
in reconstruction efficiency between non-resonant generation and a specific 
resonance channel. We also include a systematic uncertainty for the absolute 
tracking efficiency (0.2% per difference in numbers of tracks) and for K s re- 
construction efficiency (7.1% per Kg). Finally, we include a small contribution 
to the systematic uncertainty due to Monte Carlo statistics. All systematic un- 
certainties are added in quadrature. 



4 The D° -> KgKgKV Decay Mode 



The D° final state K s K s K ± 7r^ i has not been observed by previous experi- 
ments. While the D° decay modes D° -> K°K°K + ir~, D° -> K K°K-n + , 
and D° — > K°K°K + 7r~ can result in a final state K s K s K ±r K T , the decay 
mode D° — > K°K°K + ir~ is doubly Cabibbo suppressed and we assume that 
this mode is small compared to the other two Cabibbo favored modes. To 
distinguish the two Cabibbo favored decays we use the charge of the soft pion 
in the decay sequence D* + — > D°tt + . First, we demonstrate a signal for the 
two channels combined without using a soft pion tag and then we argue that 
we have evidence for independent observations in both channels. In Fig. lb 
we show the invariant mass plot for D° — > K s K s K ± tt t at L/a > 4. We fit 
the plot with a Gaussian for the signal and a 2 nd degree polynomial for the 
background using a maximum likelihood fit and we find 57 ± 10 events. 

We find no evidence for resonant substructure in D° — > K s K s K ± tc t and use a 
non-resonant Monte Carlo simulation to obtain the efficiency As a systematic 
check we find a 2.6% difference in efficiency from a D° — > a (980) + fr*(892)~ 
Monte Carlo simulation. 

Correcting for efficiency, we find: 

T(D° -> K%K%KV) 
v , L=-2 —L = 0.0106 ±0.0019 ±0.0010. 1 



We have included systematic uncertainty contributions from split samples (0), 
fit variations (0.0007), possible resonant decay (0.0001), Monte Carlo statis- 
tics (0.0001), absolute tracking efficiency (0.0001), and variation between K s 
reconstruction categories (0.00075). 

While we have established that we see a new decay channel, in order to sepa- 
rate D° -> K° s K s K + tt- and D° -> K s K s K~n + we need to tag a D° by using 
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Fig. 1. Invariant mass distribution for various D° final states: (a) Reconstructed 
mass of D° -> K s ir + ir-. There are 22835±248 events with a sigma of 11.4 MeV/c 2 . 
(b) Reconstructed mass of D° — > l£gl£gl{r = ' = 7r = '\ There are 57±10 events with a sigma 
of 5.0 MeV/c 2 . (c) Reconstructed mass of D° -» JfgiTgirg. There are 170±26 events 
with a sigma of 5.6 MeV/c 2 . 

the soft pion from the D* + — > D°tt + decay. In Fig. 2c we present the mass 
difference M^K^K^^s) - M(K s K s K ± n T ) in which we find 14.1 ±4.4 
signal events when we cut around the D° mass. Fig. 2a (2b) shows the mass dif- 
ference histogram for the events in which the soft pion has the opposite (same) 
charge as that of the kaon. The two histograms show a signal of 7.2 ± 3.4 and 
6.8 ± 2.9 events, respectively. The fits in all three histograms are performed 
by fixing the mass and width to the values returned from the Monte Carlo. 



5 The D° -> K S K S K° S Decay Mode 



Although this mode is Cabibbo allowed, it requires either final state interac- 
tions or VF-exchange to occur. Due to the limited phase space and the need 
to identify three Kg candidates, the signal can be observed without the need 
for a D* + tag or a Lja^ cut. This significantly improves the reconstruction 
efficiency. In Fig. lc we find a signal of 170 ± 26 events. 

We find no evidence for resonant substructure and use a non-resonant Monte 
Carlo simulation to obtain the efficiency A D° — > ao(980)°.K"<j simulation is 
used as a measure of the systematic variation (the difference in efficiency is 5.6 
%). For the branching ratio of this mode with respect to the D° — > K°it + it~ 
mode we find: 

r(D° — > K°K°K°) 

^— =4 S = 0-0179 ±0.0027 ±0.0026. (2) 
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Fig. 2. Mass difference between M^K^K^ir,) and M^K^K^) where n s 
is a pion from the primary vertex. Plots a, b, and c are for cases where the ir s charge 
relative to the kaon charge is opposite, same, or either, respectively. A mass cut of 
1.85 < M(KgKgK ± 7r =F 7T s ) < 1.88 GeV/c 2 to select D° events is applied. 

We have included systematic uncertainty contributions from split samples (0), 
fit variations (0.0007), possible resonant decay (0.0005), Monte Carlo statis- 
tics (0.0001), absolute tracking efficiency (0.0001), and the variation between 
K s reconstruction categories (0.0025). Several groups [5,6,7,8] have reported 
measurements of D° — > K S K S K S . 



6 The D° -> K° S K S Decay Mode 



The D° — > K°K° decay is expected to occur primarily via two VT-exchange 
diagrams, which are expected to cancel out in a four-quark model. In the 
standard six-quark model the difference between the two amplitudes is tiny, 
proportional to SU (3) flavour breaking effects, and we thus might expect the 
branching fraction for the decay to be very small (< 3 x 10~ 4 ). However, a 
Standard Model based calculation [10] predicts a relatively large branching 
fraction due to final-state interaction effects, leading to a branching ratio of 
B(D° -> K°K°) = \B(D Q -> K + K~) = 0.3%. A recent investigation [11] 
has focused on the s-channel and the t-channel one particle exchange (OPE) 
contributions. While the s-channel contribution, taken into account through 
the poorly known scalar meson / (1710), gives a small contribution, the one 
particle t-exchange gives higher contributions, with pion exchange being the 
highest. In the factorization limit [12], the branching fraction is expected to 
be equal to zero. Non-factorizable contributions in factorization-type models 
have been recently studied [13] and predict a branching fraction of about 
10~ 4 . Several experiments [14,15,16,17] have reported measurements of D° — > 
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This channel is similar to D° — > KgK§K§ in that no L/ctl cut is possible, but 
the channel requires an additional D* tag requirement. We select D* + candi- 
dates by requiring that the reconstructed mass difference AM = M D »+ — M D o 
lies within 2 MeV/c 2 of the nominal mass difference [18] of 145.42 MeV/c 2 . 

The dominant background sources for D° — > KgK^ decays are from non- 
resonant if°7r + 7r~ and -k + it~it + it~ decays. To reduce feedthrough from D° — > 
KgH + 7r~ where the 7r + 7r~ invariant mass falls in Kg mass region or from D° — > 
7r + 7r _ 7r + 7r _ where both 7r + 7r~ pairs can be misidentified as K s candidates, we 
make the additional requirement that the K s candidates reconstructed with 
silicon strip informations have a decay length significance {Ljoi) greater than 
12. We also apply a | cos^ x °l < 0-8 cu t, where 9 K o is the angle between the 
K s momentum in the D° rest frame and the D° laboratory momentum. 

In Fig. 3a the K S K S invariant mass distribution is shown for the events satis- 
fying these cuts. For the signal shape we use a double Gaussian (two Gaussians 
with the same mean) for the D° signal (79 ±17 events), as suggested by Monte 
Carlo studies, and a first order Chebyshev polynomial for the background. The 
double Gaussian shape is fixed to the one obtained from Monte Carlo simu- 
lation. In Fig. 3b the AM = D* + — D° mass difference distribution in the 
D° mass region is plotted. A Gaussian is used to fit the signal while for the 
background we use the functional form 

a(AM-m 7r ) 1/2 + 6(AM-m 7r ) 3/2 (3) 



where m n is the pion mass, the first term is from a non-relativistic model of 
phase space and the second term is the first-order relativistic correction to 
the non-relativistic model. Consistent yields are obtained from the two fits. 

The events of the normalization mode D° — > K s tt + -k' have been selected using 
the same reconstruction and fitting techniques as for D° — > K S K S mode and 
similar selection criteria to minimize systematic uncertainties. 

We measure 

T (D° -> KgK g) T(D° -> K°K°) 



T (D° -> K s n+n-) T(D° -> K°n+n- 



0.0144 ±0.0032 ±0.0016.(4) 



This accounts for the unseen i^£7r + 7r _ and K\K\ decays (note that D° — > 
K S K° L is forbidden by CP conservation). Different systematic sources have 
been investigated and the value reported includes contributions from fit vari- 
ations (0.0013) and K° s reconstruction (0.0010). 
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Yield = 79+17 



M(K S °K S « 



; 120 - 



2.030 
GeV/c ! 




Yield = 77+17 



0.136 0.146 0.156 0.166 0.176 0.186 0.196 
M(D"*-D ) GeV/c 2 




Yiekfc 113+21 



' 1.75 1.8 1.85 1.9 



M(KsKsJC 



GeV/c 



Fig. 3. Invariant mass distribution for various D° final states: (a) Recon- 
structed mass of D° —> K S K S with a D* + — D° mass difference cut. There are 
79 ± 17 events with a sigma of 12.5 MeV/c 2 . (b) Reconstructed mass difference 
of D* + — D°;D° — ► KgKg. There are 77 ± 17 events in the mass difference plot 
demonstrating consistency with (a), (c) Reconstructed mass of D° 
There are 113 ± 21 events. 



7 The D Q 



K° s K° s n 



ix Decay Mode 



This channel is Cabibbo-suppressed and was first observed by the ARGUS 
Collaboration [19]. It has a larger background than the D° — > KgKgK + ir~ 
channel due to more phase space and two pions. In addition one must eliminate 
7r + 7r~ combinations which are consistent with the Kg mass. To increase the 
signal to noise, a L/cr > 6 cut is applied. In Fig. 3c we present the KgKgir + ir~ 
invariant mass. The figure is fit with a Gaussian for the D° signal plus a 2 nd 
order polynomial. 

As the resonant structure for this channel is not known, we use a non-resonant 
Monte Carlo simulation to compute the branching fraction. 
A D° — > K* (892) ~ K* (892) + simulation is used as a measure of the systematic 
variation (the difference in efficiency is 7%). Again, we calculate the branching 
ratio of this mode with respect to the D° — > K°ii + it~ mode. 



T(D° -> K° s K° s ir + 



7T 



0.0208 ±0.0035 ±0.0021 (5) 



We have included systematic uncertainty contributions from split samples (0), 
fit variations (0.0013), possible resonant decay (0.0005), Monte Carlo statistics 
(0.0002), absolute tracking efficiency (0.0001), and the variation between Kg 
reconstruction categories (0.0015). 
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Table 1 

D° — > KgKgX Branching Fractions 



Decay Mode This Experiment PDG 2004 [18] 



0.0106 ± 0.0019 ± 0.0010 
, -, 0.0179 ± 0.0027 ± 0.0026 0.0154 ± 0.0025 

r ff°"^5"f U) 0.0144 ± 0.0032 ± 0.0016 0.0119 ± 0.0033 



s s 



r(J»°^J C|K°7r+7r ) ^ 020g ± ^ 0Q35 ± ^ Q021 ^ 31 ± ± 



8 Conclusions 



We summarize our branching ratios and compare them in Table I to the Par- 
ticle Data Group world averages where available [18]. 

We have investigated and measured the branching ratios of several decay 
modes of the charm meson D° into final states containing at least two K s 
relative to D° — > K°7t + tt~ mode. We have evidence for two new Cabibbo fa- 
vored D° decay modes, namely D° -> K° s K s K + n- and D° -> K° s K° s K-n+ 
and we have reported new measurements for D° — > K° S K S K S) D° — > K° S K S) 
and D° -> K° s K s n + n-. 
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